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Ecuaciones diferenciales de orden 
superior. Algunas aplicaciones. 


l. Sistemas masa-resorte 


1.1. Movimiento libre no amortiguado 


Consideremos un resorte de longitud | suspendido verticalmente de un soporte rígido. Si 
colgamos de él una masa m, el resorte se alargará una cantidad, llamada elongación, que 
denotamos por Al (véase Figura 1). 


Figura 1: Sistema masa-resorte. 


Las fuerzas que están actuando sobre la masa en cada momento son, por un lado, la fuerza 
de la gravedad, cuya magnitud es P = mg, donde g es la aceleración gravitatoria terrestre, y 
por otro, la fuerza restauradora del muelle que viene descrita por la Ley de Hooke: un resorte 
ejerce una fuerza restauradora Opuesta al movimiento y proporcional al desplazamiento del 
resorte; su magnitud es F, = kAl, k > 0. 

El sistema alcanza la posición de equilibrio cuando las dos fuerzas que están actuando 
son iguales; entonces se tiene que 


mg = kAl. (1) 


Si ahora desplazamos la masa una cantidad x hacia abajo desde su posición de equilibrio y 
soltamos, la fuerza restauradora será 


F, =k(x + Al). 


Considerando que los desplazamientos a partir de esta posición de equilibrio son positivos 
hacia abajo, la resultante de las fuerzas que están actuando sobre el sistema, en ausencia de 
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otras fuerzas, viene dada por 
F = mg — k(x + Al). (2) 


Por otro lado, la Segunda Ley de Newton establece que la fuerza neta aplicada sobre 
un cuerpo es proporcional a la aceleración que adquiere dicho cuerpo, siendo su masa la 
constante de proporcionalidad. Si denotamos por x(t) la posición de la masa en el instante 
de tiempo t contada a partir de su posición de equilibrio, para la que se tiene que x = 0 (ver 
Figura 1), esta fuerza viene expresada por 


F=muat). (3) 
Dado que ambas fuerzas, (2) y (3), son iguales, 
mz” (t) = mg — k(x(t) + Al) (4) 


y, teniendo en cuenta la expresión (1) se llega a la siguiente ecuación diferencial lineal de 
segundo orden: 
mz"(t) + kx(t) =0. (5) 


Esta ecuación diferencial describe el movimiento armónico simple o movimiento l:- 
bre no amortiguado. Para resolver esta ecuación de segundo orden homogénea, obtenemos 
las raíces de la ecuación característica 


mA +k=0, 


pe 
m 


puesto que k > 0. La solución general de esta ecuación diferencial vendrá dada por 


x(t) = e, cos (12) + ca sen (Es) (6) 


que, escrita en términos de la frecuencia (número de oscilaciones por unidad de tiempo) 


k 
w=w/—, queda escrita como 
m 


x(t) = c, cos(wt) + ca sen(wt). (7) 


que son 


El movimiento descrito por la expresión (7) es un movimiento oscilatorio periódico, con 
: ; 27 m pe : 
amplitud constante y cuyo periodo es T' = — = 271,/—. El movimiento del sistema masa- 
W k 
resorte es una vibración de frecuencia constante que no depende de las condiciones iniciales, 


tan solo depende de k y m. 


Ejemplo 1.1 Si situamos una masa de 5 kg en un resorte, éste se alarga 10 cm. Liberamos 
la masa 8 cm por debajo de la posición de equilibrio. ¿Cuál es la ecuación del movimiento 
suponiendo un movimiento armónico simple? (Tómese el valor aproximado de y = 10 m/s?). 
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De la expresión (1) obtenemos que 


mg  5-:1000 


OE 


= 500 kg - s ?. 


El P.V.I. que describe la dinámica del resorte, para este caso, vendrá dado por la ecuación 
diferencial (5) junto con las condiciones iniciales dadas, es decir, 


Las soluciones de la ecuación característica, 54? + 500 = 0, son 
A= TO, 
y por tanto, la ecuación que describe el movimento es 
x(t) = c, cos(10t) + cz sen(10t). (8) 
Imponiendo las condiciones iniciales en 


cy cos(10t) + cz sen(10t), 
—10 c, sen(10t) + 10 cz cos(10t), 


OS 

A 

aa 
po] 


tenemos que 
8. =-2(0). => 6N 
0 == 2(0)-=: 1063, 


y, por tanto, las constantes c, y cz toman los valores c; = 8 y cz = 0. La solución particular 
única para este ejemplo viene dada por 


x(t) = 8 cos(10t). 


En la Figura 2 se muestra la gráfica de esta función x(t), en la que se aprecia el movimiento 


217 
periódico, con periodo 1' = 10 y amplitud 8. 


87 


Figura 2: Movimiento libre no amortiguado. 
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1.2. Movimiento libre amortiguado 


El caso del movimiento libre no amortiguado no es muy realista. Si suponemos, por 
ejemplo, que el resorte está sumergido en un fluido (aire, agua, ...), actuará sobre él una 
fuerza amortiguadora que es proporcional a la velocidad del cuerpo y con sentido opuesto; 
esto es, 


F. =-—cx(t), 


donde c > 0 es una constante conocida como coeficiente de amortiguamiento. Así, en la 
expresión (4) habría que considerar también este término, quedando ahora escrita como: 


mz" (t) = mg — k(x + Al) — cx'(t). (9) 


Procediendo como en el caso anterior llegamos a la siguiente ecuación diferencial que describe, 
en este caso, el denominado movimiento libre amortiguado: 


mx"(t) + cx'(t) + kx(t) =0. (10) 


Se trata de una ecuación diferencial lineal de segundo orden homogénea. La ecuación 
característica se escribe como 


mA +cA+k=0 (11) 


cuyas raíces son 
—c + ve? — dkm 


2m 


A= 


Dependiendo del signo de c? — 4km distinguiremos tres casos: 


(a) Si c? — 4km > 0, tendremos dos raíces reales distintas, 


, —=c+ ye? — 4km ) =c=— ye — km 
Li y IS A 
2m 2m 


La solución de la ecuación diferencial (10) viene dada por 
at= ae tae tz, (12) 


Puesto que los valores A; y A2 son negativos, Jim x(t) = 0, es decir, la posición de la masa 
—00 


tiende a su posición de equilibrio. Además, se prueba fácilmente que la masa puede pasar 
a lo sumo una vez por su posición de equilibrio. Se dice que este sistema está sobreamor- 
tiguado. Esto ocurre cuando el coeficiente de amortiguamiento, c, es grande comparado con 
la constante del resorte, k. 


Ejemplo 1.2 Supongamos que una masa de 1 kg alarga 5 m un resorte. Determinar la 
ecuación del movimiento libre amortiguado si la masa se libera dos metros por debajo de 
la posición de equilibrio con una velocidad ascendente de 3 m/s suponiendo que la fuerza 
amortiguadora es 3 veces la velocidad instantánea. 
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De la expresión (1) obtenemos que 
mg 1-10 155 
k= ===>——=2kg:s “. 
Aa cd 
Dado que c= 3, el P.V./ para este ejemplo de movimiento libre amortiguado es 
x"(t) + 3x2'(t) + 2x(t) =0, 
2 0=32, 
x'(0) = —3. 
Las soluciones de la ecuación característica, A? +31 +2 = 0, son 
Ar =-—1 y Ao = -—2 


y por tanto, la ecuación que describe el movimento viene dada por 


x(t) = cet + ee”. (13) 
Imponiendo las condiciones iniciales sobre 
x(t) = cqet+oe”, 
xt) = —ejet—2ce*, 


se llega al siguiente sistema de ecuaciones para C, y Ca: 
2 = x(0) = c1+a>, 
33 = x(0) = —c,-—2cz, 


La solución particular que describe el movimiento del 


cuya solución es cy = 1 y ca = l. 


sistema viene dada por 
x(t) =et+e*, 


En la Figura 3 se muestra la gráfica de x(t) en la que se aprecia como la masa tiende a su 


posición de equilibrio. 


Figura 3: Sistema sobreamortiguado. 


(b) Si c? —4km < 0, tendremos como soluciones de la ecuación característica (11) dos raíces 


complejas conjugadas, 
—c + iv4km — c? —c— iv4km — e? 


Aj = A = 
d 2m y ñ 2m 
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La solución de la ecuación diferencial (10) viene dada ahora por 


Lee (E ) 
(¡0 2m"C0S t 


2m 


x(t) 


, —3mt 
FC3€ 2m*" sen 


vVikm — ec? ) 


2m 


e (E ) (E )) 
e 2m" | (1 COS t] + casen | ————t ; 


2m 2m 


El movimiento descrito por x(t) es un movimiento oscilatorio que se va amortiguando con 
el tiempo ya que su amplitud, que depende de e72*, tiende a cero cuando t > 00. Cuando 
ocurre esto, se dice que el sistema está subamortiguado, y se produce cuando el coeficiente 
de amortiguamiento, c, es pequeño comparado con la constante del resorte, k. 


Ejemplo 1.3 Determinar la ecuación del movimiento de un sistema masa-resorte para el 
caso m =1kg, c=2N s/m y k= 10N/m suponiendo que la masa se libera desde la posición 
de equilibrio con una velocidad descendente de 3 m/s. 


El P.V.T que tenemos que resolver este ejemplo, de acuerdo con la expresión (10) y las 
condiciones iniciales del enunciado, viene dado por 


x"(t) + 2x'(t) + 10x(t) =0, 
z(0)=0, 
2 (0=3: 
Su ecuación característica asociada, A2 +2, + 10 = 0, tiene por raíces 
Su solución general será, por tanto, 
x(t) = cie cos(3t) + cae” sen(3t) = e** (c, cos(3t) + ca sen(3t)). 
Las condiciones iniciales impuestas en 


e”* (e, cos(3t) + ca sen(3t)) , 
—e* (c, cos(3t) + cz sen(3t)) + e? (—3c, sen(3t) + 3c2 cos(3t)) , 


NS 
8.8 
== 
pol 


permiten llegar al siguiente sistema de ecuaciones para C¡ y Ca: 


O. = "m0 <= cr 
3 = (0) = —c +3c2, 


de donde c; =0 y c2 = 1. La solución particular que describe el movimiento del sistema se 
escribe, entonces, como 
xa(t) = e*sen(3%). 


La Figura 4 muestra la gráfica de x(t). Nótese que el movimiento es oscilatorio en torno a 
la posición de equilibrio con amplitudes cada vez menores que tienden a cero. Cabe destacar 
que este comportamiento no depende de las condiciones iniciales impuestas. 
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Figura 4: Sistema subamortiguado. 


(e) Si c? — 4km = 0, tendremos como solución de la ecuación característica (11) una raíz 


real doble, A = e En este caso, la solución de (10) viene dada por 
m 


x(t) = cre + cate? = e% (e, + cat). (14) 


El valor de Á es negativo, y por tanto, ón x(t) = 0. De la expresión (14) se deduce, además, 
=00 


que la masa no podría pasar más de una vez por su posición de equilibrio. Este sistema se 
dice que está críticamente amortiguado ya que una pequeña variación en la fuerza de 
amortiguamiento o en la masa haría que el sistema estuviese sobreamortiguado o subamor- 
tiguado. 


Ejemplo 1.4 Una masa de 2 kg se sujeta a un resorte cuya constante es k = 2N/m. 
Supongamos que sobre el sistema está actuando una fuerza amortiguadora que es igual a 4 
veces la velocidad instantánea. Determinar la ecuación del movimiento si la masa se libera 
1 m por debajo de la posición de equilibrio con una velocidad descendente de 1 m/s. 


El P.V.I que describe el movimiento del sistema masa-resorte en este caso será 


=1 


21" (t) + 41 (1) + 2x(t) =0, 
J=1 
£.(0)= L 


Su ecuación característica, 21? + 44 +2 = 0, tiene por solución A= —1 doble, y por tanto, 
la solución general de la ecuación diferencial vendrá dada por 


x(t) = ce *+cxte*=e* (c, + tez). 


Imponiendo las condiciones iniciales en 


et (c, +tc2), 
e * (01 +tco) +e cz, 


OS 

po 

a 
Al 
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llegamos al siguiente sistema de ecuaciones para C, y Ca: 
LES ON 
1. =. (0) 


de donde se obtiene que cy = 1 y cz =2. La solución particular que describe el movimiento 
del sistema se escribe, entonces, como 


C1, 
=C1 + Ca, 


x(t) =e*(1+2t). 


En la Figura 5 se representa la gráfica de x(t). 


Figura 5: Sistema críticamente amortiguado. 


Ejemplo 1.5 Consideremos el mismo sistema masa-resorte del ejemplo anterior con las 
condiciones iniciales dadas por x(0) =1 y x'(0) = —3. 


En este caso, la solución particular viene dada por la expresión 
x(t) = e *(1—2t). 


Podemos observar que en instante t = 1/2, la masa pasa por su posición de equilibrio, y sólo 
en ese instante (véase Figura 6). 


Figura 6: Sistema críticamente amortiguado. 
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1.3. Movimiento forzado amortiguado 


Consideremos ahora que una fuerza externa, F(t), está actuando sobre el sistema. Si en 
la expresión (9) añadimos el término correspondiente a esta fuerza, tendríamos que 


mu" (t) = mg — k(x+ Al) —cx (0) + F(0, (15) 


y, de modo análogo a los casos anteriores, llegaríamos a que la ecuación diferencial que 
describe el movimiento del sistema, que en este caso se denomina movimiento forzado, 
queda escrita como 

max" (t) + cx'(t) + kx(t) = F(t). (16) 


Así, tenemos una ecuación diferencial lineal de segundo orden no homogénea. Estudiaremos 
como caso particular fuerzas externas periódicas de la forma 
F(t) = Fosen(£2t), o bien, F(t) = Fo cos(£2t), ya que en estos casos se producen fenómenos 
físicos interesantes. 


Ejemplo 1.6 4 un sistema masa-resorte amortiguado cuyos parámetros son m = 1kg, 
c=4N5s/m yk = 3 N/m se le aplica una fuerza externa dada por F(t) = 5 cost. Determinar 
la ecuación que describe el movimiento del sistema suponiendo que r(0) =0 y 2'(0) =0. 


El P.V.I. que resuelve este ejemplo es 


x"(t) + de (t) + 3x(t) = 5 cost, 
(0) =:0, 
PA 07=0: 


Primero resolvemos la ecuación diferencial homogénea asociada 2%, (t) +42 (6) +324(t) =0. 
Su ecuación característica 
AM+4A+3=0 


tiene como soluciones A; = —1 y A2 = —3, por lo que la solución general de esta ecuación 
viene dada por 

ult) = ce * + ce *. 
Para encontrar una solución particular de la ecuación diferencial completa aplicaremos el 
método de los coeficientes indeterminados. Ensayamos, entonces, con una solución particular 
de la forma 


xp(t) = Acost+ Bsent 


de modo que 
Tp(t) =—Asent + Bceost, 
xp(t) =—A cost — Bsent. 


Sustituimos ahora las expresiones de <p(t), vp(t) y vp(t) en la expresión de la ecuación 
diferencial completa y tenemos 


—Acost — Bsent + 4(—Asent + Bcost) + 3(Acost + Bsent) = 5cost. 
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Operando e igualando los coeficientes de sent y cost llegamos al siguiente sistema de ecua- 
ciones en Á y B: 
2A+4B = 5, 
l 4A+2B = 0, 


1 
cuya solución es A== y B=1. Así, 


1 
uplt) = a cost + sent 
y, por tanto, la solución general de la ecuación diferencial completa se escribe como 
—( —-3t 1 
x(t) = ce *+ce "+ 3 cos t + sent. 


A continuación, para determinar los valores de c, y cz imponemos las condiciones iniciales 
en 


1 
x(t) = cet +ee + 3 cost + sent, 


x'(t) = —cye”* — 3cge* — , sent + os t, 
con lo que tenemos 
1 
0 = x(0) = Coca + 
= (0) = —c -—3c2+1, 
j 5 3 de j . el 
de donde se obtiene que c, = qe Ca = de La solución particular que describe el movimiento 
del sistema es E A ] 
x(t) = Se q + 7 Cost + sent. 


Figura 7: Movimiento forzado amortiguado. 


Nótese que, según avanza el tiempo, únicamente se mantiene el movimiento oscilatorio 
generado por la aplicación de la fuerza externa (ver Figura 7). 
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1.4. Movimiento forzado no amortiguado 


Este caso particular del movimiento forzado tiene especial interés ya que puede dar lugar 
a un fenómeno conocido como resonancia y que describiremos más adelante. 

Puesto que estamos suponiendo ausencia de amortiguación, c = 0, la expresión (16) se 
escribe ahora como 


mx" (t) + kx(t) = F(t). (17) 


Supongamos, como se comentó en el caso anterior, que la fuerza externa es periódica con 
una frecuencia £2 y viene dada por la expresión F(t) = Fo cos($2t), donde Fy € R. Por tanto, 
tendremos 


mz" (t) + ka(t) = Fo cos($9t). (18) 


La ecuación diferencial homogénea asociada, mx, (t) + kug(t) = 0, corresponde a la que 
describe el caso del movimiento armónico simple cuya solución escrita en términos de la 


frecuencia, w = 4 /—, venía dada por la expresión (7): 
m 


xp(t) =c, cos(wt) + ca sen(wt). 


Para resolver la ecuación diferencial completa, aplicando el método de los coeficientes 
indeterminados, ensayaremos con la solución particular 


xp(t) = Acos(12t) + Bsen(S2t), 


de modo que 


p(t) = —ASN sen(12t) + BN cos(12t), 
x(t) = —A? cos(12t) — B2? sen(£t). 


Sustituyendo <p(t) y sus derivadas en la expresión (18) se tiene 
m (—AS2? cos($2t) — BS2 sen(£2t)) + k(Acos($2t) + B sen(S2t)) = Fo cos(S2t). 


Operando e igualando los coeficientes de cos((£2t) y sen(£2t) llegamos al sistema de ecuaciones 
en Ay B 


B(-m12 + k) =0, 
Fo 


de donde se obtienen los valores A = Ema Y B =0. Por tanto, el movimiento del sistema 
—-m 


( A(=m12 + k) = Fo, 


viene descrito por 


x(t) = c, cos(wt) + ca sen(wt) + 7 cos(12t). (19) 


Eo 
k— m£ 
Para interpretar este movimiento resulta conveniente expresar la amplitud de las oscila- 


ciones producidas por la aplicación de la fuerza externa, +———5, en función de la frecuencia 


k — m422” 
de las oscilaciones del movimiento libre, w = 4/—, esto es 
m 


x(t) = c, cos(wt) + casen(wt) y a cos((2t). (20) 


12 Ecuaciones Diferenciales Ordinarias de Orden Superior 


Nótese que este análisis sólo es válido si w 4 S2. Sin embargo, de (20) se puede observar que 
cuando £2 —> w la amplitud de las oscilaciones tiende a infinito, lo que da lugar al fenómeno 
conocido como resonancia (ver Figura 9). 


Ejemplo 1.7 Al sujetar una masa de 2 kilogramos a un resorte cuya constante es 
k = 32 N/m, éste queda en reposo en la posición de equilibrio. A partir de t = 0, se aplica al 
sistema una fuerza externa dada por F(t) = 4cos(2t). Encontrar la ecuación del movimiento 
en ausencia de amortiguación. 


El P.V.I. para este caso es 


2x"(t) + 32x(t) = 4cos(2t), 
x(0) =0, (21) 
x'(0) =0 


Como siempre, resolvemos primero la ecuación homogénea asociada, 2x%,(t) + 3224 (t) =0. 
Para ello buscamos las raíces de la ecuación característica 24? + 32 = 0 que son A¡ = 4i y 
A2 = —41, por lo que su solución general se escribe como 


xp (t) = cc, cos(4t) + cz sen(4t). 
Suponemos ahora que una solución particular es de la forma 

xzp(t) = Acos(2t) + B sen(2t), 
de modo que 


xp(t) = —2A sen(2t) + 2B cos(2t), 
ap (t) = —4A cos(2t) — 4B sen(2t). 


Sustituyendo las expresiones de xp(t) y v(t) en la ecuación diferencial completa (21), 


2(-4A cos(2t) — 4B sen(2t)) + 32(A cos(2t) + B sen(2t)) = 4 cos(2t), 


dl 
operando e identificando coeficientes, obtenemos que Á = 6 y B =0. Por tanto, la solución 


general de la ecuación completa vendrá dada por 
1 
x(t) =c, cos(4t) + ca sen(4t) + A cos(2t). 
Imponemos ahora las condiciones iniciales en 


1 
x(t) = c, cos(4t) + ca sen(4t) + 6 cos(2t), 


1 
x'(t) = —4c, sen(4t) + 4cz cos(4t) — 3 sen(2t), 


esto es, 
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Figura 8: Movimiento forzado no amortiguado. 


1 
de donde se obtiene que c; = 6 y C2 =0. La solución particular que describe el movimiento 
del sistema es 1 y 

a cos(4t) + 6 cos(2t), 


cuya representación gráfica se muestra en la Figura 8. 


Ejemplo 1.8 Consideremos el mismo sistema masa-resorte del ejemplo anterior al que se 
le aplica ahora una fuerza externa dada por F(t) = 3 cos(4t). 


El P.V.I. que hay que resolver en este caso es 


2x"(t) + 32x(t) = 3 cos(4t), 
£(0)=0, (22) 


para el que se tenía que la solución general de la ecuación homogénea asociada venía dada 
por 
xy(t) = cc, cos(4t) + cz sen(4t). 


Puesto que, en este caso, la fuerza externa, F(t) = 3 cos(4t), es una función que es solución 
particular de la ecuación diferencial homogénea, para encontrar una solución particular de 
la ecuación completa ensayaremos con 


xp(t) = At cos(4t) + Btsen(4t), 
de modo que 


Yp(t) = Acos(4t) — 44t sen(4t) + Bsen(4t) + 4Bt cos(4t), 
xp (t) = —8A sen(4t) — 164t cos(4t) + 8B cos(4t) — 16Btsen(4t). 


Sustituimos ahora estas funciones en la expresión (22), 
2(-84A sen(4t) —164t cos(4t)+8B cos(4t) -16Bt sen(4t))+32(4t cos(4t)+Bt sen(4t)) = 3 cos(4t), 


y procediendo de manera análoga a los ejemplos anteriores obtenemos los valores A = 0 y 


BR 16" de modo que la solución general de la ecuación (22) resulta 


3 
x(t) = c, cos(4t) + ca sen(4t) + 16 sen(4t). 
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Imponemos las condiciones iniciales en 


3 
x(t) = c1 cos(4t) + ca sen(4t) + 16% sen(4t), 


3 3 
x'(t) = —4c, sen(4t) + 4cz cos(4t) + 16 sen(4t) + q cos(4t), 


y obtenemos los valores c, = c2 = 0. La solución particular que describe el movimiento del 
sistema viene dada por 


ai) a sen(4+). 


Como puede observarse en la Figura 9, en este caso se produce el fenómeno de resonancia, 
ya que la frecuencia de la fuerza externa coincide con la del movimiento oscilatorio libre. 


Figura 9: Movimiento forzado no amortiguado. Resonancia. 


Ejercicios propuestos 


1. Supongamos que una masa de 2 kg alarga 5 m un resorte. Determinar la ecuación del 
movimiento libre amortiguado si la masa se libera 2 m por encima de su posición de equi- 
librio sin velocidad inicial suponiendo que la fuerza amortiguadora es 4 veces la velocidad 
instantánea. 


Solución: x(t) = —2e*(cost + sen t) 


2. Consideremos una masa de 5 kg sujeta a un resorte de constante k = 20 N/m. Sobre este 
sistema está actuando una fuerza amortiguadora de constante c = 20 N s/m. Si la masa se 
suelta 2 m por debajo de su posición de equilibrio con una velocidad inicial descendente de 
1 m/s, determínese la ecuación del movimiento del sistema. 


Solución: x(t) = 2e% + 5te=* 


3. Determinar la ecuación del movimiento de un sistema masa-resorte formado por una masa 
de 8 kg suspendida de un muelle de constante k = 32 N/m suponiendo que, en ausencia de 
amortiguación, actúa una fuerza externa dada por F(t) = 16 cos(4t). La masa se libera desde 
la posición de equilibrio con una velocidad inicial ascendente de 1 m/s. 
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1 1 1 
Solución: x(t) = 6 cos(2t) — a sen(2t) — 6 cos(4t) 


4. A un sistema masa-resorte no amortiguado cuyos parámetros son m = 1 kg y k = 9 N/m se 
le aplica una fuerza externa dada por F(t) = 6sen(3t). Determinar la ecuación que describe 
el movimiento del sistema suponiendo que x(0) =1 y x'(0) =0. 


1 
Solución: x(t) = cos(3t) + 3 sen(3t) — t cos(3t) 


2. Ejercicio resuelto con Maple 


Una masa de 4 kg cuelga de un muelle de constante k = 36 N/m. Una fuerza externa dada por 
3l 
F(t) =3c0s (5) está actuando sobre el sistema. Determinar la ecuación del movimiento 


del mismo suponiendo que la masa se libera desde la posición de equilibrio, sin velocidad 
inicial y que no está actuando ninguna fuerza amortiguadora. 


> restart: with(plots): 
La ecuación que modela el sistema masa-resorte es 

> ecuacion:=4*diff(x(t),t,t)+36*x(t)=3*c08 (31/10xt) ; 

ecuacion := 4 des (t) +36 x (t) = 3 cos (5 ) 
dt? 10 

Primero, resolvemos la ecuación homogénea asociada, 414,(t) + 36 xp (t) = 0. Para ello, 
obtenemos las raíces de su polinomio característico, 41? + 36 = 0. 

> solve(14*172+36=07,11)); 

[1 =31), Xi =-—31) 

Por tanto, la solución general de la ecuación homogénea tiene como expresión 

>  xH:=t->c1x*cos(3*t)+c2*sin(3xt); 

1H := tr cl cos(3t) + c2 sin (3t) 

Para encontrar una solución particular de la ecuación completa ensayamos con 


31 PO E 
xP(t) == Acos (5+) + Bsin (5) 


> xP:=t->Ax*cos(31/10*xt)+Bxsin(31/10xt); 


31 (31 
xP := to ÁAcos (50) + Bsin (50) 


Sustituimos zp(t) en la ecuación diferencial completa 
> 4xdiff (xP(t),t,t)+36x*xP(t)=3xsin(31/10)x*t; 
61 dl 61 al Sl 
—=— A —t]-—Bsinl| — =3 sin| — 
25 cos (5 ) 35 so (5 ) 3 sin (5): 
y obtenemos las constantes A y B igualando coeficientes 


> solve(1-(61/25)*A=0,-(61/25)*B=37F,1A,B); 
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75 
A paa B o 
e a) 


La solución general de la completa viene dada por 
> xg:=t->clxcos(3xt)+c2*sin(3*t)-75/61*c08(31/10*t); 


61 10 
Determinamos, ahora, el valor de las constantes cl y c2 imponiendo las condidiones iniciales; 
antes tenemos que obtener la expresión de 2'(t) 


> AR 


715 31 
rg := te cl cos (31) + c2 sin(3t) — — cos (51) 


y 465 . (31 
=3 cl sin(3t) +3c2 cos(3t) + 199 (5 ) 


>  dxg:=t->-3xc1*sin(3*t)+3*c2*c0s (3*t)+(465/122)*sin((31/10)x*t); 


465 31 
dig := t+>H-3c1 sin(3t) +3c2 cos (3t) + 193 sin (5 ) 


> solve(txg(0)=0,dxg(0)=0+,tfc1,c2)); 
75 
cl =—,c2 =0 
u=a 
La expresión de la solución particular que describe el movimiento del sistema masa-resorte 
es: 
> x:=t->75/61x*xc08(3*t)-75/61*c08(31/10*t); 


E A E aL, 
IT := 61 COS 61 COS 10 


> plot(x(t),t=0..40*Pi,color=black,thickness=1); 


ATA 
A 


Figura 10: Fenómeno de pulsación. 


o 


Este fenómeno, conocido como pulsación, (véase Figura 10) se produce cuando la frecuencia 
de las oscilaciones libres y la de las de la fuerza externa son casi iguales, esto es, w = £2, 
(véase expresión (20)). La función 


se puede expresar como 


(+) = 150 en. ÚN , 
alt) = 7 sen opt) sen la? )> 
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(tras aplicar equivalencias trigonométricas). Se observa como la amplitud de este movimiento 


, z . ; 61 ' 
oscilatorio, que viene dado por el término r5N sen (5). depende del tiempo de manera 


periódica, puesto que la solución posee dos frecuencias no proporcionales. 


3. Deflexión de una viga 


En la construcción de gran parte de las estructuras se utilizan vigas. Estas vigas sufren 
deformaciones que pueden ser debidas a su propio peso o a la acción de fuerzas externas. 
El cálculo de los valores que representan estas deformaciones es de gran importancia ya sea 
para determinarlos en puntos específicos a lo largo del eje de la viga o para comprobar que 
no exceden los límites permitidos. 

Supongamos una viga horizontal, homogénea, cuya sección transversal es uniforme en 
toda su longitud. En estas condiciones, su eje de simetría (curva que une los centroides de 
las secciones transversales de la viga en ausencia de cargas, incluida la de su propio peso), es 
una recta. La curva que describe la deformación que experimenta la viga (su desviación re- 
specto al eje de simetría) al aplicar una carga sobre ella se conoce como curva de deflexión 
o curva elástica (ver Figura 11). 


eje de simetría curva de deflexión 


Figura 11: Eje de simetría. Curva de deflexión. 


La expresión analítica de esta curva se podrá determinar, como veremos, a partir de una 
ecuación diferencial lineal de cuarto orden. En este caso, habrá que resolver un problema 
de contorno o problema con valores en la frontera, esto es, la solución particular 
se obtendrá imponiendo condiciones de la función o alguna de sus derivadas en puntos 
diferentes. Estas condiciones de contorno dependerán de cómo esté apoyada la viga, según 
se indicará más adelante. 

Consideremos una viga homogénea de longitud L a la que se le aplica una carga en un 
plano vertical en el que está contenido el eje de simetría (ver Figura 11). Supongamos que 
el eje OX lleva la dirección del eje de simetría y sea u(x) la función que describe la curva de 
deflexión, es decir, la función que en cada punto x proporciona los valores de la deflexión, 
que supondremos positiva hacia abajo. En estas condiciones, vamos a obtener la ecuación 
diferencial lineal que permitirá obtener la expresión analítica de u(x). 
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En Teoría de Elasticidad se muestra que la relación entre el momento de flexión o 
momento flector en cada punto x, M(x), y la carga por unidad de longitud, w(x), 
viene dada por 

M"(2) = —wv(x). (23) 


A su vez, se tiene que, M(x) es proporcional a la curvatura de la curva elástica, k(x), esto 
es, 


M(x) = —Elk(o), (24) 


donde E es el módulo de Young del material de la viga e / es el momento de inercia de 
una sección transversal de la viga respecto del eje neutro (ver Figura 12). Al producto El 
se le conoce como rigidez a flexión de la viga. Ambos parámetros, E e 1, los suponemos 
constantes. 


Eje neutro 


Superficie neutra 


Figura 12: Eje neutro. 


A partir de estas dos expresiones, (23) y (24), se tiene que 
wr=.EL (e: (25) 
Por otro lado, se sabe que la expresión de la curvatura de una función u(x) es 


We) 


AD = yr 


Puesto que los valores de u(x) y u'(x) son muy pequeños en las estructuras más comunes, 
podemos despreciar u/(x)? y, por tanto, aproximar la función curvatura por 


k(1) = u' (1), (26) 
de modo que la expresión (25) quedaría escrita como 


Elu" (2) =w(w), (27) 


siendo ésta la ecuación diferencial lineal de cuarto orden cuya resolución proporciona la 
expresión analítica de la curva de deflexión de la viga. 

Como se indicó anteriormente, las condiciones de contorno van a depender de cómo estén 
físicamente los extremos de la viga. Puede ocurrir que la viga esté empotrada en ambos 
extremos, en voladizo (empotrada en un extremo y libre en el otro), apoyada en ambos 
extremos o empotrada en un extremo y apoyada en otro (ver Figura 13). 
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(a) Viga empotrada (b) Viga en voladizo 


(c) Viga apoyada (d) Viga empotrada en un 
extremo y apoyada en el otro 


== | A 


Figura 13: Apoyos en vigas. 


(a) Extremo empotrado 


Cuando uno de los extremos de la viga está empotrado, se verifica que en él no hay 


desplazamiento, es decir, 
u(x) =0, 


y además, la curva de deflexión en ese punto es tangente al eje OX, es decir, 
a (1) =05 


Suponiendo que el extremo de la viga que está empotrado corresponde al valor x = xp, en 
este caso, las condiciones de contorno para ese punto serían 


l u(xo) =0, 


wire) ='0. 


(b) Extremo apoyado 


Cuando un extremo está apoyado, se verifica que en ese extremo, de manera análoga al 
caso anterior, 


la) =0: 


Además, en este punto el momento flector es cero, esto es, M(x) = 0. Teniendo en cuenta la 
expresión (24) esto es equivalente a que k(x) = 0 y según la expresión (26) podemos decir 


que 
u'(x) =0. 


Por tanto, en este caso, las condiciones de contorno en el extremo apoyado x = xy sON 


Í u(zo) =0, 


ww (Ly).=01 


(c) Extremo libre 


En el extremo libre, el momento flector es cero; es decir, 
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y el esfuerzo cortante, V(x) = M'(x), también es cero, con lo que 


uv" (2) =0. 
Por tanto, en este caso, las condiciones de contorno en el extremo libre x = xy son 


u"(£0) =0, 


(to) =U: 


Ejemplo 3.1 Calcular la deflexión de una viga de longitud L empotrada en ambos extremos 
si una carga wey está distribuida de modo uniforme a lo largo de su longitud. Determinar 


también la ley del momento flector y el esfuerzo cortante. 


En este caso se tiene que 


w(t)= ww, Va € [0, £]. 


El problema de contorno que hay que resolver es 


Elu'Y (2) = wo, 


u(0) =u'(0) =0, 
u(L) =w(L) =0. 
W 
La ecuación diferencial escrita como uV (2) = 2 se resuelve de manera inmediata inte- 


erándola cuatro veces; esto es, 


e) 7ro+ A, 
>) Wwy 1? 
3 A 2 
aula E | > EBI+C, 


4 3 2 
Wa YT Ax Bzx 
= | | LCx+D 
ul) = 54 +7 iS 


Las constantes A, B, C' y D se obtienen imponiendo las condiciones de contorno dadas 


anteriormente: 


u(0)=0 $ D=0, 


v(0)=0 $ C=0, 


woL?* ABS. Be 


L = T =0, 

AS 
aL? AL? 

n= i EBL=0, 

a(Ei=0: SEI a 0 
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de donde se obtiene que 
waL m0 wL? 
2E1' —1EP” 


Por tanto, la curva de deflexión viene definida por la función 


C=0,  D=0. 


Es 2 


Wo 4 Wwol 3 pe 
MEL - 


AE” —12EI” * 


u(z) 


T T T 
0,25 0,50 0,75 1 
x 


0,001 


Figura 14: Curva de deflexión para w=1, El =2 y L=1. 


El momento flector, M(x) = —Elu”(x) vendrá dado por 


wo a weLl pl? 
Mai== | E 
(2) O E E 


L 
Los valores M(0) = M(L) = Ma 


de gran importancia en muchos métodos de análisis de estructuras. 


se denominan momentos de empotramiento y son 


-0,083 


coc! 


Figura 15: Momento flector para w =1, El =2 y L=1. 


wal 
Y por último, el esfuerzo cortante será V(x) = M'(x) = —wyt + E 


Figura 16: Esfuerzo cortante para w =1, El =2 y L=1. 


Ejercicios propuestos 


1. Calcular la deflexión de una viga de 1 metro de longitud empotrada en el extremo izquierdo 
y libre en el derecho, si una carga “y está distribuida de modo uniforme a lo largo de su 
longitud. 
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A a z 
24EI 6EI 4El 


2. Calcular la deflexión de una viga de 1 metro de longitud apoyada en los dos extremos, si 
una carga we está distribuida de modo uniforme a lo largo de su longitud. 


Solución: u(z1) 


Ea 0 at — 2 a? y se z 
2EI 1281” ' 24EI 
3. Calcular la deflexión de una viga de 1 metro de longitud empotrada en el extremo izquierdo 


y apoyada en el derecho, si una carga wy está distribuida de modo uniforme a lo largo de su 
longitud. 


Solución: u(x) 


= wo Gi Wo ya 
24 EI 48El 


wo 


2 
16E1" 


Solución: u(zx) 


4. Ejercicio resuelto con Maple 


Hallar la expresión analítica de la curva de deflerión de una viga de longitud L empotrada en 


—zx 
el extremo izquierdo y libre en el derecho suponiendo que w(x) = we . Para los valores 


we =3El y L= 1/2 representar gráficamente la curva y calcular la defleción máxima. 


En primer lugar, cargamos los paquetes necesarios: 


> restart:with(plots): 


La ecuación diferencial que hay que resolver es u!Y (1) = 2. la definimos como 
>  u4:=x->(w0x(L-x)/(Lx*E1)); 
wO (L—x) 


Integramos cuatro veces, añadiendo en cada una de ellas las constantes de integración 
A, B,C y D correspondientes. 
1 
wO | Lx— zx? 


LEI 
> 13:=x->w0x*x(Lx*x-(1/2)*x72)/(L*EID)+A; 


1 
wO ES ¿a) 
EA 


LEI 


> intlul1(x),x); 


A 


> intluB3(G),x); 
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> 12:=x->w0x((1/2)*L*x"2-(1/6)*x73)/(L*ElI)+Ax*x+B; 


UY i=-> 


> intlu2GD),x); 


TA 


> intílul(x),x); 


il 
dr zx +5 Ba" +Cz 


UN= ho 


1 
| +2 Ba? y ED 
pe y BT Czx 


Una vez obtenida la solución general de la ecuación diferencial, u(x), aplicamos las condi- 
ciones de contorno correspondientes al caso de una viga empotrada en un extremo y libre en 
el otro: 


> solve(fu(0)=0,u1(0)=0,u2(L)=0,u3(L)=07$,1A,B,C,DF); 
en 1 w0L _1w0 E? 
NN A e 


.0=0,D=0) 


> subs(%,u(x)); 


0 df 4 ll 5 
A o ET 1 w0Lx3 1 w0L?x? 


LEI 12 Ad, 2 - El 


> up:=x->w0x((1/24)*L*x"4-(1/120)*x75)/(Lx*E1) -(1/12)*w0x*L+*x"3/El+ 
(1/12) *w0x*L72x*x"2/ElI; 
1 
Ol —La*- —a? 
AS (5 id 0”) 1 w0La3 1 w01?%x? 


pie LEI =D En “Y El 


La expresión analítica de la función que describe la curva de deflexión viene dada por: 
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1 1 
EA do A 
ES (5 “19 ) O O 
LEI 12 El E 1 


I 
en la que ahora sustituiremos los valores de “y y L dados en el enunciado, 


up(1) = 


> subs((w0=3*E1,L=1/2,up()); 
A O A E 


> flex:=x->(1/8)*rx"4-(1/20)*x7"5-(1/8)*xx73+(1/16)*x72; 


1 1 1 1 
fer i=ae creta eta q? 
8 20 8 16 


y obtenemos para este caso concreto la función 


1 
flen(a) => ¿au 1 E lr 


Su representación gráfica es: 
> plot(-flex(x),x=0..1/2,color=black) ; 


Obsérvese que dibujamos la función flex(x) con signo negativo puesto que se considera el 
eje OY positivo hacia abajo. 


El valor de la deflexión máxima se alcanzará en el extremo de la viga; esto es: 


> evalf(flex(1/2)); 
0.006250000000 


Ejercicios propuestos 


1. Hallar la expresión analítica de la curva de deflexión de una viga de longitud L apoyada 
en ambos extremos suponiendo que w(x) = wyz?/L para los valores Wy = 4E1 y L = 2. 
Representar gráficamente la curva. 


Solución: u(1) = A A A 
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2. Hallar la expresión analítica de la curva de deflexión de una viga de 1 metro de longitud em- 
potrada en su extremo izquierdo y apoyada en el derecho suponiendo que 
w(x)=(2+4cosx) El. 


1 93 49 
Solución: u(1) = od +4cos 1 + (s cos 1 — a) a+ (5 — 7 cos 1) 7-4. 


3. Hallar la expresión analítica de la curva de deflexión de una viga de longitud L empo- 
trada en el extremo izquierdo y libre en el derecho suponiendo que w(x) = wesen(rx/L) para 
los valores y = 3El y L = 1/2. Representar gráficamente la curva y calcular la deflexión 
máxima. 


3 sen(211) 1 3% 32 
16  m 47 87 81 
Deflexión máxima: 0.01384720518. 


Solución: u(1) = 
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